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Abstract 

The short-time evolution of a growing interface is studied within the framework of 
the dynamic renormalization group approach for the Kadar-Parisi-Zhang (KPZ) equation 
and for an idealized continuum model of molecular beam epitaxy (MBE). The scaling 
behavior of response and correlation functions is reminiscent of the "initial slip" behav- 
ior found in purely dissipative critical relaxation (model A) and critical relaxation with 
conserved order parameter (model B), respectively. Unlike model A the initial slip ex- 
ponent for the KPZ equation can be expressed by the dynamical exponent z. In 1+1 
dimensions, for which z is known exactly, the analytical theory for the KPZ equation is 
confirmed by a Monte-Carlo simulation of a simple ballistic deposition model. In 2+1 
dimensions z is estimated from the short-time evolution of the correlation function. 

PACS numbers: 68.35 Fx, 64.60 Ht, 05.70 Ln, 05.40 +j 
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I. Introduction 



Interface formation and growth are typical processes in nonequilibrium systems. From a 
technological point of view two important examples are fluid flow in porous media (oil in 
rock) Q and deposition of atoms during molecular beam epitaxy (MBE) Q |2|. It is ex- 
pected that at times much later than typical aggregation times and on macroscopic length 
scales these interfaces develop a characteristic scaling behavior, where the scaling expo- 
nents fall into certain dynamic universality classes jl], ^, [| (see below). In certain cases, 
however, interfaces can also show turbulent, i.e., spatial multiscaling behavior [Q]. Usually 
a d-dimensional interface is embedded in d + 1-dimensional space such that the interface 
position at time t can be described by a height function h(x.,t), where x denotes the lat- 
eral position in a d-dimensional reference plane given by, e.g., the surface of a substrate 
in MBE. Complete information about the scaling behavior is contained in the dynamic 
structure factor, which is related to the time displaced height-height correlation function 
C(x — x.',t,t') = (/i(x, t)h(x! , t')) — (/i(x, t)) (/i(x', t')), where a laterally translational invari- 
ant system is assumed. For t, if — ► oo and finite \t — if\ the correlation function displays the 
asymptotic scaling behavior 

C(x - x', t, = |x - x'\ 2a F c (\t - i'|/|x - x'H, (1-1) 

where a denotes the roughness exponent and z is the dynamic exponent [l], §]• F° r a laterally 
translational invariant system the interfacial width w 2 (t) = (h 2 (x,t)) — (/i(x, t)) 2 is only 
a function of t and displays the scaling behavior w(t) ~ t@ for late times, where [3 = 
a/z is the growth exponent. For MBE as an example the scaling behavior displayed in 
Eq. (|l.l| ) gives access to the exponents a and z both experimentally by reflection high energy 
electron diffraction (RHEED) (see, e.g., chaper 16 of Ref.|j|]) and theoretically by continuum 
models Q and Monte-Carlo simulations |^, g|. Since the advent of the scanning tunneling 
microscope (STM) direct imaging techniques for interfaces have also become an important 
experimental tool 0. 

Continuum descriptions of interfacial growth processes can be obtained from general 
symmetry principles and conservation laws obeyed by the growth process Q. The resulting 
coarse grained growth model is given by an evolution equation for h(x, t) which has the form 
of a Langevin equation with Gaussian distributed noise. This has first been done in Ref.|7]] 
for the sedimentation of granular material and leads to the well known Edwards- Wilkinson 
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(EW) equation. It is given by 

|/t(x,t) =z.V 2 / l (x,t)+r ? (x,t), (1.2) 

where the noise 7y(x, t) has a Gaussian distribution with (ry(x, t)) = and 

(r/(x, t)r/(x', 0> = 2D<5(x - x')£(* - *')• (1-3) 

The parameters i/ and -D are assumed to be constants and averages (. . .} are taken over the 
noise distribution. From explicit solutions of Eq. (|1.2| ) the exponents z and a are known 
exactly in any dimension d of the interface: 

z = 2 and a = (2-d)/2. (1.4) 

Note that d = 2 is the critical dimension of Eq.(|1.2j). One has a < for ci > 2 so that the 



height-height correlation function C(x — x', £, if) according to Eq.(l.l) decays with increasing 
distance |x — x'|. In d = 2 (a = 0) the correlations increase logarithmically. 

The simplest possible nonlinear extension to the EW equation was first considered sys- 
tematically in Ref.[§|. The resulting Langevin equation, usually denoted as the Kadar-Parisi- 
Zhang (KPZ) equation, is given by 

f h(x, t) = i/V 2 /i(x, t) + |(V/i(x, t)) 2 + r/(x, t) (1.5) 

with Gaussian distributed noise according to Eq.( |1.3[ ). The additional parameter A is again 
assumed to be a constant. In the long time limit Eq.( |1.5[ ) has a global symmetry which is 
commonly denoted as Galileian invariance. This invariance originates from the equivalence 
of Eq.(|1.5|) to the Burgers equation for a vorticity-free velocity field v(x, t) = — V/i(x, t) and 
can be stated as follows. If h(x,t) solves Eq.Ql.5D for some noise function rj(x,t) then 



h'(x, t) = h(x - wt, t) - iw • x + ^w 2 t (1.6) 

is a solution of Eq.(|l.5|) for the noise function r/'(x, t) = r/(x — wt, t) and any constant vector 
w. An important consequence is that the exponents z and a of the KPZ equation fulfill the 
exact scaling relation [||, |l0| 

a + z = 2. (1.7) 

Note that in d = 2 the EW exponents also obey Eq.( |l.7[) . The exponents of the KPZ equation 
are exactly known only in d = 1, where 

z = 3/2 and a = 1/2 (1.8) 
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due to the existence of a dissipation fluctuation theorem |], [Tl]]. In d = 2 numerical inves- 
tigations indicate z oi 1.6 and a ~ 0.4 [0. For d > 2 the asymptotic scaling behavior is 
either governed by the EW exponents (see Eq.( |1.4[ ), weak coupling regime) or by another 
set of exponents inaccessible by analytical methods (strong coupling regime) depending on 
the value of the effective coupling constant g = DA 2 /(4z/ 3 ) [|], [9|, [lOtl - In d = 3 numerical 
evidence suggests z ~ 1.7 and a ~ 0.3 in the strong coupling regime fil still indicating rough 



interfaces in contrast to EW scaling behavior in d = 3 (see Eq.(1.4)). Furthermore, it is 
interesting to note that the nonlinearity in Eq.( |1.5| ) is the most relevant one, i.e., if present it 
renders all other nonlinearities irrelevant in the renormalization group sense in the long-time 
limit. For intermediate times, however, the presence of other nonlinearities in the Langevin 
equation gives rise to various crossover phenomena [jl], [l2|]. The EW equation and the KPZ 
equation for A / thus represent two different universality classes for interfacial growth. 



For A < Eq.(1.5) can be viewed as a model for interface corrosion rather than growth ||. 

With special regard to MBE growth it is worth noting that the requirement of mass 
conservation in ideal MBE |l^] explicitly excludes the KPZ nonlinearity from a corresponding 
coarse grained continuum theory. A simple Langevin equation for ideal MBE has been 



proposed in Ref.O] (see also Refs.|14, |15|]): 



|/t(x, t) = -^V 4 /i(x, t) + Ai V 2 (V/i(x, t)) 2 + ??(x, t), (1.9) 



where rj(x, t) is chosen according to Eq. fll.3p . Mass conservation in combination with Eq.(1.3) 
immediately leads to the exact scaling relation 2a — z + d = for Eq. (|1.9| ). Furthermore, a 
global symmetry analogous to Eq. (|1.6| ), which can be written in the operator form [|l(| 



x^x-2wtV 2 , /i^/i-^-w-x (1.10) 



for any infinitesimal vector w, yields the second exact scaling relation a + z = 4 fl3| , 16]. The 
exponents z and a for ideal MBE are therefore known exactly in any dimension of physical 
interest: 

z = (8 + d)/3 and a = (4-d)/3 (1.11) 



indicating d = 4 as the critical dimension of Eq.( |1.9| ). 

In this paper Eqs. (|l.5| ) and ( |l.S| ) are used as paradigms for continuum descriptions of 
interfacial growth processes. In linear theory (i.e. A = Ai = 0) their dynamical exponents are 
given by z = 2 and z = 4 (see Eqs.( |l.4[) and ( A.13j )), respectively, and therefore Eqs. (|1.5[) and 
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(|1.9| ) may be viewed as nonequilibrium analogues of the dynamical models A and B for critical 
relaxation, respectively. In order to invesitgate the scaling behavior of, e.g., C(x — x',t,t') 
for t' <C t the initial condition h(x, t = 0) = motivated by deposition processes is used 
simultaneously with Eqs.( |1.5| ) and fll.3| ) or Eqs.(1.9) and ( |1.3|) , respectively. Perturbative 
and nonperturbative aspects of short-time scaling for the two models are discussed in Sees. II 



and III within the framework of dynamic renormalization jlT], 15, 19]. Numerical results 
from ballistic deposition are presented in Sec. IV and a summary of the main results is given 
in Sec.V. 



II. KPZ equation 

Due to the spatial translational invariance of the deposition processes studied here calcu- 
lations are most conveniently performed in Fourier space. With the definition h(x, t) = 
(27r)~ rf J d d q exp(iq • x)/i(q, t) for the Fourier transform the dynamic functional J~[h,h] for 
the KPZ equation [||, [l?], [2(| can be written as the sum of the Gaussian part 

Jo[h, h] = J J™ dt |r>^(q, £)M-q, t) - h(q, t) (^("Q, *) + WK~<\, *)) } t 2 - 1 ) 

and the interaction part 

J 1 [h,h]=-^J ^ J ^ J^dtqx-qzhi-qx-qarfhiqirfhfat), (2.2) 



where h(q, t) is the Fourier transform of the response field IS ] . The initial condition /i(q, 0) 



0, which is implicitly assumed in Eqs,([2.1|) and ( |2 . 2| ) , breaks the temporal translational 
invariance of the KPZ dynamics. In a more general form this broken symmetry can be 
expressed in terms of an additional contribution to J7b which is localized at the time "surface" 
t = 0: 



Uh] = I l 7^b(Mq,0) -/* (q)) 2 . (2.3) 



2 7 ( -0(M q ,O)-/ io (q)) 2 . 
From the analogy of Eq. (|2.3| ) with surface contributions to the Ginzburg Landau functional 



in the theory of static surface critical phenomena [21] and dimensional arguments the only 
possible fixed point values of c under the renormalzation group are c = ±oo and c = 0. 
In the latter case additive renormalizations of c are supposed to be absorbed in c itself 
which can be accomplished by the dimensional regularization scheme. On the other hand 
Eq.( |2.3| ) generates a distribution function exp(—J' s [h]) of initial configurations h(q, 0) of 
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the deposition field which leaves the fixed point value c = oo as the only choice due to the 
requirement of normalizability of distribution functions. Deviations of c from this fixed point 
value therefore only generate corrections to scaling |2l| which will be disregarded here. 
From Eq.( |2.3|) one then has the initial condition /i(q, 0) = /io(q)- As shown in Appendix A 
/io(q) can be incorporated into a source contribution to the dynamic functional (see Eqs.( [A3| ) 
and ( |A.6[ )) and therefore we stick to /io(q) = in the following. The correlation and the 



response propagator are now easily derived from Eq.(2.1). The results are summarized in 
Appendix A. 

The introduction of an initial condition striktly speaking also breaks Galileian invariance 
(see Eq.(|1.6|)). If one demands /i(x, 0) = as the initial condition for h then /i'(x, t) solves 
Eq.( |1.5| ) with the new initial condition h'(x,0) = — w • x/A. However, as indicated above 
one only has to transform the source fields accordingly in order to restore the old initial 
condition. Therefore the Galilei transformation (see Eqs.( |1.6|) , ( |A.3 ), and ( |A.6| )) 



h'(q, t) = e-^Hq, t) - (2vr)^w • ^«5(q) , fc'(q, t) = e^'Mq, t) , 
j'(q, t) = e~^ 1 j(q, t) , j'(q, t) = e"^ j(q, t) + (2vr)^w • ^<5(q)<5(i) 



(2.4) 



restores the Galileian invariance of the generating functional so that the corresponding Ward 
identities (see RefiQ) remain valid. Note that Eq. (|2.4] ) should be read as an infinitesimal 
transformation, i.e., terms of order w 2 have been neglected. 

The renormalization group treatment of Eq.([E^) can now be set up following standard 
procedures ||, |9|, ffifl . For the case at hand it is most convenient to combine the dimen- 
sional regularization scheme for the KPZ equation B with the treatment of the short-time 



singularites documented in Ref.[17|. One defines the effective coupling constant 

g = D\ 2 /(A^) (2.5) 

and the renormalized parameters u R , D R , and u ||, |2^] 

v R ee Z v v , D R ee Z D D , and u = Z g gfi £ /(2 d " 1 vr a!/2 (2 - d/2)), (2.6) 

where e = d — 2 and fi is an arbitrary momentum scale which absorbs the naive dimension 
of g (see Eq. fl2.5| )). One finds the renomalization factors |^, (see also Appendix C) 

Z u = l + ^ + 0(u 2 ), Z D = l-^ + 0(u 2 ), 

Z h = Z h = 1 , and Z g = Z D Z~ 3 , (2.7) 
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where the 1/g poles indicate the presence of ultraviolet singularities ||. The nonrenormal- 
ization of h and h indicated in Eq.( |2.7D is exact and a consequence of Eq.( B.l[ ) (see Appendix 
B). The relation Z g = Z]jZ~^, which is equivalent to X R = A, is a consequence of Galileian 
invariance in the long-time limit (see Eq.(1.6) and Refs.[||, |9|, ^0|) and therefore also holds 
to all orders in perturbation theory. The renormalization group flow at late times is then 
governed by the Wilson functions |9|, ^] 

( U ( U ) = ^ U + 0( U 2 ), ( D ( u ) = - u + 0(u 2 ), 

/3(u) = (d-2 + $ D (u)-3(„(u))u, (2.8) 

where the relation between f3(u), Cd(u), and ( u (u) is again exact. The higher order correc- 



tions to and Cd indicated in Eq^^) vanish in d = 1 due to the existence of a fluctuation- 
dissipation theorem || The fluctuation-dissipation theorem also requires Z u — Zy) in 
d = 1 so that Cu(u) = Cd(u) and v/D = v R /D R (see Eqs.(|l|) and tfufy). We also want 
to emphasize here that ( u (u) and (d(u) as given by Eq. (|2.8|) , like any other finite order 



perturbation theory, do not give access to the strong-coupling regime of Eq.( |1.5|) for d > 2. 

In analogy with critical phenomena in semiinfinite geometries |2l| modifications of the 
scaling behavior of response and correlation functions must be expected in the "time surface" 



t = [17]. In order to determine the corresponding anomalous short-time scaling dimensions 
of response and correlation functions we introduce two new renormalization factors Zq and 
Zq by the renormalization prescription (see also Ref.|TT 



/i(q, 0) = zl /2 h R {(\, 0) and h(q, 0) = zj /2 ^(q, 0). (2.9) 



These Z-factors are determined by Eq^f^l) and the operator identity 

^h(ci,t = 0) = 2Dh(ti,t = 0) (2.10) 
derived in Appendix B. For the weak coupling regime of the KPZ equation (d = 1) the 



perturbative analysis of Appendix B consitutes a rigorous proof of Eq. Q2.10 ) and the relations 
which follow from it (see below). In the strong coupling regime (d > 2), however, the 
corresponding perturbative analysis does no longer provide a rigorous proof of Eq.( |2.10| ), 
because relations which are valid order by order in perturbation theory may be violated at a 
strong coupling fixed point (see Appendix B). This has to be kept in mind for the following 
considerations, although the perturbative result can be regarded as evidence in favor of the 
general validity of Eq,( |2.lti| ). 



7 



From Eq.(B.l) for if = we immediately find the exact identity Zq = 1. Insertion of 
Eqs.( |2.6| ) and ( |2.9| ) into Eq.(2.10) leads to the second exact identity Zq = Z^ 2 Zq which 
determines Zq in terms of the known Z- factor Zp (see Eq.(2.7)). These identities translate 
into the exact relations 

Co(«) = and Co(n) = -2C D (n) (2.11) 

among the corresponding Wilson functions (see also Eq.( p.8[ )). From Eq,( 2.1l| ) one concludes 
that (i) the response function G(q, t, t') does not exhibit an anomalous scaling dimension in 
the short-time limit t' — » (i.e., t' <C t) and that (ii) the anomalous short-time exponent 
of the correlation function Co(q,i, t') can be expressed by long-time exponents (see Eq.( |1.8j ) 
and the following text). These properties set KPZ short-time dynamics markedly apart from 
model A. 

In order to determine the short-time scaling exponent of C(q, t, t' <C t) we employ the 
"short distance expansion" |l7]] /i(q, t' — > 0) = a(t')-^j /i(q, t' = 0) + . . . inside the correlation 
function C which means that 







C(q, t, t' « t) = o{t')— C(q, t, t' = 0) + . . . . 
Employing the renormalization prescriptions given by Eqs.([2.6|) and ( |2.9| ) one finds 



(2.12) 



a(t') = Z 1/2 a R (fi,t',u) and ^ C(q, t, t' = 0) = Z%*j^ C H (fi, q, t, t' = 0, «) (2.13) 

for the corresponding renormalized short distance expansion (see Eq.( 2.12j )). Using dimen- 
sional analysis the renormalized functions defined by Eq,( |2.13| ) can be written in the scaling 
form 







A/2 



d 



d 

— C R {n, q, t, t' = 0, u) = D(n)g(x, y, u) 



with 



with 



y' = v{n)li t' and 



(2.14) 



q/^ and y = ^(/u)/i 2 t, 



where /i has been chosen as the renormalization group flow parameter. It is now straightfor- 
ward to derive the renormalization group equations for the dimensionless scaling functions 
f(y',u) and g(x,y,u) defined by Eq.( 2.14| ). Using Eqs. (|2.6|) and ( |2.9| ) one obtains 

Co(«) 



(2 + Uu))y'-^+(3(u) ° 







-x- + (2 + an)y- +/3 (n)|- + C^) + 



du 2 
Co(u) 



f(y',u) 
g{x,y,u) 



and 



(2.15) 



At the infrared stable renormalization group fixed point u = u* Eq.( |2,15| ) has the solutions 

f(y',u*)=y'^ an d g(x,y,u*)=y- (2vD+71o)/2z g'(x z y), (2.16) 

where r] a = Ca(u*) for a = u, D, 0, z = 2 + rj v , and g' is a scaling function left undetermined 
by Eq,( p,15l ). Combining Eqs.( 2TT^ ), ( |2.14 ), and ( |2,16 ) one finds after a few manipulations 



C(q, t, t' « t) = {t'/t) l+V £ hr- z fc(h\ z t) (2.17) 

for the short-time scaling behavior of the correlation function. For u* = one obtains the 
EW scaling exponents (see Eq. Ql.4D ) and tjq = in Eq.( |2,17| ). For any nonzero fixed point 
u* the exact scaling relation rj£> = 3z — 4 — d holds (see Eq.( |2.8p ), which is equivalent to 
Eq. Ql.7p . From Eq.( p.li ) one finally obtains for the short-time exponent (see Eq.( 2.17] )) 



1 + 77o/(2z) =0 = 1- ri D /z = (d + 1)1 z - 2. (2.18) 

In d = 1 the exact value 9 = 4/3 can be obtained from Eq.(|l. 8;). From numerical estimates 
for z in d = 2 and d = 3 (see Sec. I) one obtains 9 ~ 1.7 and 9 ~ 2.1, respectively. The 



exponent relation given by Eq.( 2.18 ) simply means that the short-time and the long-time 
scaling behavior of the correlation function are identical, i.e., the short-time scaling behavior 
can be obtained by extrapolating the t'-dependence of C(q, t,t') from t' ~ t to t' = 0. 



In fact, the scaling relation given by Eq. Q2.18 ) can be derived independently by analyzing 



the fluctuation spectrum of the interface displacement velocity averaged over a macroscopic 



portion of the interfacial area 22]. 

Finally, we remark that some alternative scaling forms for C can be obtained from the 
definition of the growth exponent (3 = a/z which leads to 9 = d/z + 2(5. The scaling 
behavior displayed in Eq,( |2.17j ) can then be written in the simplified form C(q, t,t' <C t) = 
t' d 9c{\<l\ z t), where gc{y) = U~ G fc{y)- I n real space the correlation function has the scaling 
form C(x,t,f < t) = (t'/t) e \x\ 2a G c (t/\x\ z ). 

The absence of anomalous scaling exponents for G(q, t, t') for t' <t does not neccessarily 
mean that G is analytic for t' — > 0. Exponents describing the asymptotic short-time behavior 
are in general functions of the dimensionality d and therefore may take noninteger values for 
certain d. Similar considerations apply to the crossover behavior of G for t — > oo with fixed 
t — t'. For details we refer to Appendix C, where some results from perturbation theory are 
discussed in the case d = 1. 
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III. Ideal MBE 



In terms the deposition field /i(q, t) and the response field /i(q, t) the dynamic functional 
J[h,h] for Eq.flD]) @, 0, Q is also written as the sum of the Gaussian part 

Jo[h, h}= J dt |^(q, t)fc(-q, t) - h(q, t) (J^K-q, t) + ^(q 2 ) 2 /i(-q, t)) } 

(3.1) 

and the interaction part 
J 1 [h,h}=\ 1 J — -p J dt(qi +q 2 ) 2 qi -q 2 /i(-qi -q 2 ,t)/i(qi,t)/i(q2,0> (3-2) 



where the initial condition fo(q, 0) = is again implicitly assumed in Eqs.( |3.1| ) and (3.2). As 
described in Sec. II and Appendix A this special initial condition is sufficient to study the 
short-time scaling behavior of response and correlation functions for ideal MBE. The results 
of Gaussian theory as implied by Eq.(|3.1|) are summarized in Appendix A. 



The further analysis of Eq.(1.9) can be carried out along the lines of the analysis of the 
KPZ equation presented in Sec. II. First, we note that the invariance under the infinitesimal 
transformation given by Eq.( l.lCj ) in presence of the initial condition /i(q, 0) = is restored 
by the transformation 



h'(q, t) = e 2 "^ wt Mq, t) - (2^w • A^q) , fc'( q , t ) = e 2 ^ wt h(q, t) , 
j'(q, t) = e 2 ^^ w *j(q, t) , j'(q, t) = e 2 ^^j(q, t) + (2vr) d ^-w • A<S(q)*(*)> 



(3.3) 



where terms of the order w 2 have been neglected. In analogy with the Galileian invariance 
of Eq.( |1.5| ) this symmetry leads to the nonrenormalization of the nonlinearity: X R = Ai 



(see also Eq.( [2.7D and Refs.|13|, ^6|). Second, Eq.(1.9) has the global symmetry of mass 
conservation which in contrast to Eq. ( |2.6| ) leads to the additional nonrenormalization of the 
noise correlation amplitude (see Eq.( |1.3| )): D R = D [13, 16]. If one defines an effective 
coupling constant by |L3| 

gx = DXj/uf (3.4) 
and the renormalized parameters v R , D R , and u 

u^Z^, D« S Z D D, and u S ™ , (3.5) 
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where e = d — 4 and /i is an arbitrary momentum scale which absorbs the naive dimension 
of g\ (see Eq.( |3.4D ) then the renomalization group results for Eq.( |l,9| ) in the long-time limit 
can be summarized as follows (see also Appendix C): 

Zux = 1 + + 0(u 2 ) , Z D = 1, 
Z h = Z h = 1 , and Z g = Z^ 3 . (3.6) 

The 1/e poles indicate the presence of ultraviolet singularities and the nonrenormalization 
of h and h indicated in Eq. Q3.6| ) is again a consequence of Eq.( |B.l| ). The corresponding 
renormalization group flow is therefore governed by only two nontrivial Wilson functions, 
namely 

( Vl (u) = ^u + 0(u 2 ) and 0{u) = (d - 4 - 3^ (u))u, (3.7) 

where the relation between (5{u) and C,vi{ u ) is exact (see Eq. (|3.6|) ). For any infrared stable 
fixed point u* ^ Eq.(j37^) yields Q Vl (u*) = rj Vl = (d — 4)/3 from which the exponents given 
by Eq. fll.li ) follow directly. 



In order to investigate the short-time behavior of the response and the correlation function 
of Eq.( |1.9| ) short-time renormalization factors Zq and Zq are defined as in Eq.(p.Sj). From 
Eqs.( |B.i] ) and (2.10), which also hold for Eq.( |1.9|) (see Appendix B), one immediately obtains 



Zq = 1 and Zq = Z^ 2 Zq = 1, where Eq.(|Q|) has been used. We thus conclude that 
in contrast to KPZ dynamics for ideal MBE neither the response function G(q,t,t') nor 
the correlation function C(q, t, t') exhibit anomalous scaling behavior for if <C t which is 
reminiscent of the short-time behavior of model B in critical relaxation |p.7||. Finally we note 



that in contrast to model B the noise in Eq.( |1.9[ ) is not conserved (see Eq.(1.3)). Eq. Ql.9| ) 
with purely conserved noise has been first considered in Ref.(l(| (see also Ref.|j]]). The 
qualitative short-time behavior is the same as described here. However, with special regard 
to MBE the case of purely conserved noise does not play the same central role as Eq.( |1.9| ) 
with nonconserved noise Q{ and we therefore refrain from discussing any details here. 

Concerning the asymptotic short time behavior of G(c[,t,t') and C(q, t, t') and the 
crossover to their asymptotic long-time behavior one finds properties which are similar to 
the KPZ behavior mentioned in Sec. II. Some details obtained from perturbation theory are 
reported in Appendix C. 
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IV. Ballistic deposition 



The scaling behavior of C(q, t,t' <C t) according to Eq.( 2.17| ) can be tested numerically by 



a Monte-Carlo simulation of a simple ballistic deposition model on a lattice with periodic 
boundary conditions . For convenience we restrict ourselves to d = 1 here. The continuum 
description used in Sees. II and III is replaced by a discretized description according to 

h(x,t) = h(x = aj,t = n/(FL)) = ahj(n), (4.1) 

where the lattice constant a is assumed to be the same both in the plane of the substrate 
and perpendicular to it. The lattice has L sites, F is the incoming particle flux, and n 
is the number of deposited particles. Furthermore, the incoming particle flux F has been 
normalized to unity, so that t in Eq.( |4.lD is dimensionless and given by the number of 
deposited layers. Finally, hj(n) defined by Eq.(|4.1j) is also dimensionless and denotes the 
number of particles deposited at lattice site j after n particles have been deposited on the 
lattice. Ballistic deposition on a one-dimensional substrate is defined by the deterministic 
growth rule 

hj(n + 1) = max(/jj_i(n), hj(n) + 1, hj+i(n)) (4-2) 

(see, e.g., RefiQ), where the site j in Eq.( [4.2[ ) has been selected randomly from the L sites 
of the lattice. For periodic boundary conditions hi(n) and h^n) are treated as nearest 
neighbors in Eq.( |4.2| ). 

In order to measure the scaling behavior of C(q, t,t') given by Eq.( [2,17 ) for the above 
discrete model a discrete Fourier transform is defined by 

1 L 2-7T 

h q {n) = j j Y, hj(n)e- iqaj with q = —m, (4.3) 



where m is an integer and < m < L — 1. Using Eq.(4.3) we define the discrete version of 



the height-height correlation function in Fourier space by 

C L (?,t 1 t / ) = ((^(n)-<^(n)>)(^(n / )-(^(n / )>)), n = FLt, and n' = FLt', 

(4.4) 

where the triangular brackets (...) denote an average over different realizations of the de- 
position process and the time arguments t and t' are reintroduced for convenience. For 
the measurement of the short-time exponent 6 (see Eq.( |2,l§| )) it is sufficient to measure 
CL(q,t,t') for q = 0. In this case Eq.(4.4) defines the time displaced correlation function of 
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Fig. 1: Correlation function Cl(0, t, t') in d = 1 as a function of t' jt for 0.002 < 
t'/t < 1 and L = 480 (dotted line), L = 960 (dashed line), and L = 1920 (dash- 
dotted line). The errorbars are shown only at a few selected points in time and 
represent one standard deviation. The solid line displays a power law with the 
theoretical short-time exponent 9 = 4/3. The data follow this power law rather 
accurately in the interval 0.03 < t'/t < 0.4 (see main text). 



the spatially averaged deposition height h q= o(n), which can be measured very quickly during 
the simulation. In practice a measurement is done after the deposition of one layer, i.e., the 
time step is At = 1. 

Like a real deposition process the simulation is characterized by an a priori unknown 
microscopic aggregation time t a . A scaling behavior of Cl according to Eq, (|2.17j ) can only 
be observed for t' 3> t a . On the other hand t' <C t is required for Eq.( [2.17| ) to hold, so that 
short-time scaling is restricted to the time window t a <C t' <C t. Furthermore, the lattice 
size L must be chosen sufficiently large in order to avoid the onset of finite-size crossover 
effects if i' 1 / 2 ~ L when t' is still much smaller than t. For the simulation described here 
t = 2000 and L > 480 fulfill the above requirements. In order to cope with the very small 
signal to noise ratio in each measurement of Cl(0, t,t') for t' <S t averages are taken over 
10 5 realizations. These are distributed over 40 individual runs at every point in time for 
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all lattice sizes. The result is displayed in Fig.[l], where Cx(0, t,t') is shown as a function of 
t' jt for fixed t = 2000 and for L = 480, 960, and 1920. For clearity the statistical error is 
shown only at a few points in time. As can be seen from Fig.|l| there is slightly more than 
one decade in t' jt available to determine the short-time exponent 9. Using the least square 
method in the interval 0.03 <t'/t< 0.4 one finds 

9 = 1.349 ± 0.005 (4.5) 

for L = 1920 as the best estimate for 9 from the data shown in Fig.Q. Although the agreement 
with the theoretical value 9 = 4/3 is very good, there is still a systematic deviation well 



outside the statistical error, which is one standard deviation in Eq.(4.5). One source of 
systematic errors is the finite lattice size. For example one finds 9 = 1.37 for L = 480 and 
for L = 240 (not shown in Fig.[l]) one even has 9 = 1.40 which indicates that finite lattice 



corrections to Eq.( |2.17 ) are still visible in Eq.(4.5) as a small systematic deviation of 9 from 



its theoretical value. Furthermore, Eq.(2.17) displays only the leading singular behavior of 
the correlation function in the KPZ universality class. For the ballistic deposition model 
studied here corrections to scaling not captured by Eq.( 2.17| ) may lead to sizeable numerical 



deviations. Therefore, the exponent 9 measured here should be interpreted as an effective 
exponent. However, the numerical data for Cl(0, t, t') follow a simple power law governed by 
this effective exponent quite accurately. Deviations from this power law begin to show only 
for t' jt > 0.4, where one is clearly outside the short-time limit, and for t' jt < 0.03, where 
microscopic aggregation effects come into play. 

In d = 2 the ballistic deposition model described here can be used to estimate the 
dynamic exponent z of the KPZ equation. The growth rule for ballistic deposition on a two 



dimensional square lattice with L x L lattice sites is the natural extension of Eq.(4.2): 



hj,k(n + 1) = max(hj-i ik (n), hj )k ~i(n), hj ik (n) + 1, h j+1;k (n), h jjk+1 (n)), (4.6) 

where periodic boundary conditions have been assumed. The lattice momentum has two 
components and is given by q = (27r/L)(mi, 7712), where m-i and mi are integers and < 
mi,?7i2 < L — 1. The correlation function Cx(q, t,t') is defined as in Eq.(4.4), where the 



lattice Fourier transform /i q (n) of the deposition field is defined in analogy with Eq.( |4.3j ). 
Note that n = FL 2 t with F normalized to unity relates n and t in this case so that t is again 
given by the number of layers deposited on the substrate. The short-time exponent 9 can 
be measured as described above by measuring Cx(q = 0,t, t') (see Eq.(|4.4D) for t' t. In 
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order to keep the amount of CPU time needed for the simulation within reasonable limits we 
reduce t to t = 1000 and take averages over 2 x 10 4 realizations of the deposition process. It 
turns out that a linear lattice size of L = 120 sites is already sufficient to uniquely identify 
at least one decade for the scaling variable t'/t in which Ci(0,t,i') obeys the simple power 
law predicted by Eq.( 2.17| ). The overall behavior of CL(0,t,t') for L > 120 is qualitatively 



the same as displayed in Fig.|T] so that we refrain from reproducing it here. For L = 240 and 
0.01 < t'/t < 0.1 we obtain 

9 = 1.655 ± 0.052 (4.7) 

from a least square fit as the best estimate for from the available data. Using Eq.( [2,18| ) 
we obtain the estimate 

z = 1.642 ± 0.052 (4.8) 



from Eq.(4.7) as our estimate for the dynamical exponent z in the KPZ universality class in 
d = 2. A corresponding estimate for z can be obtained for L = 120 which differs by less than 
half a standard deviation from the value given by Eq.(4.8), so that finite size effects can be 



neglected within the statistical error. Finally, we note that according to Eqs.( |1.7| ) and 
one has a = 0.358 ± 0.052 for the roughness exponent. These values are in agreement with 
other numerical data for z and a in d = 2 (see chapter 8 of Ref. [0 for a collection of recent 
estimates) and they therefore provide some support for the general validity of Eqs,( 2?LC| ) and 

V. Summary and discussion 

The following main results have been obtained: 

1. The short-time dynamics of the KPZ equation can be analyzed in close analogy to 
the short-time behavior of model A in critical relaxation. Starting from the operator 
identity given by Eq.( [2.10] ) the analogy can be summarized schematically: 

f /i(q, 0) = 2Dh{^ 0) => Z Z = Zl 2 Z Z 

model A: KPZ: (5.1) 

Z D = (Z/Z) 1 / 2 Z = Z = 1 

Zq = Z Q ^ 1 Z = 1 , Z = Z D 2 . 
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In contrast to model A the anomalous short-time scaling dimension 9 of the deposition 
field h is given by the dynamical exponent z (see Eq.( [2,18|) ), whereas the response field 
h does not exhibit an anomalous short-time scaling dimension. The capability of ana- 
lytical methods with regard to a full quantitative description of the crossover behavior 
from short to long times for the KPZ equation is limited. A perturbative analysis 
combined with dimensional considerations indicate that q 2 (i — t') 2 /t d / 2 is the scaling 
argument which governs the leading finite-time corrections to the asymptotic long time 
scaling behavior of the response function G(q, t, t') in d = 1. In the correlation function 
C(q = 0,t, t') finite-time corrections persist indefinitely. A quantitative description of 
the full scaling behavior can probably be obtained by combining perturbative methods 
with mode coupling theory |1C]. 



2. The short-time dynamics of ideal MBE according to Eq.( |1.9| ) can be analyzed in close 
analogy to the short-time behavior of model B in critical relaxation. Starting again 
from Eq.( |2.10| ) neither the deposition field h nor the response field h exhibit anomalous 
short-time scaling dimensions, which is the same behavior as observed for model B [[17|]. 
In contrast to the KPZ equation the infrared stable renormalization group fixed point is 
finite in any dimension of physical interest. Therefore purely perturbative methods can 
be used to investigate the short-time to long-time crossover behavior of the response 
and the correlation function within, e.g., an e = d — 4 expansion. In combination with 
dimensional arguments the perturbative analysis indicates that the scaling argument 
q 4 (t — t') 2 jt d l^ governs the leading finite-time correction to the asymptotic long-time 
behavior of the response function G(q, t, t'). In the correlation function C(q = 0, t, if) 
finite-time corrections again persist indefinitely. 



3. With a simple ballistic deposition model Eq. (|2.18|) can be used to measure the dy- 
namical exponent z for the KPZ universality class from a simulation of the short-time 
behavior of the height - height correlation function. Although such a simulation in 
principle requires short computer times, the overall benefit is somewhat limited due to 
the small signal to noise ratio in the correlations for t' <C t (see Fig.|l]) which in turn 
must be compensated for by running the simulation with high statistics. In d = 1, 
where Eq.( |1.8[ ) gives the exact scaling exponents, the numerical results for 8 agree very 
well with the theoretical value 9 = 4/3 which is equivalent to z = 3/2. In d = 2 
Eq,( [2.18|) has been successfully used to obtain a numerical estimate for the dynamical 



16 



exponent z in the KPZ universality class (see Eq.(4.8)). 

Finally, it should me mentioned that the short-time scaling behavior of the magnetization 
in an Ising model with model A (Glauber) dynamics can be efficiently used to determine 



the dynamic and static critical exponents in the Ising universality class [ 23 [ . It would be 
interesting to see, to what extent Monte-Carlo methods similar to those described here and 
in Ref.|23j can be used to study the asymptotic long-time scaling behavior of interfacial 
growth models from their short-time dynamics. The scaling relation between 6 and z may 
also open an alternative path for direct numerical investigations of the KPZ equation. 
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Appendix A: Gaussian theory 

The Gaussian part J7o of the dynamic functional for Eq. (|1.5| ) is the same as for model A of 
critical relaxation [fn]] and can be written in the symmetric form 

= \lwfC dt (M-^)M-*,i)) A J ^ J I , (A.l) 

where /i(q, 0) = and the response field fulfills the additional condition /i(q, oo) = 0. The 
self- adjoint matrix operator A is then given by 

/ 2D -4 -z^q 2 \ 

A = \ m . (A.2) 

In terms of the source fields j and j introduced by adding the source term 

Jjlkh] = J 4SdJ™ dt (Mq,t)J(-q,t) + Mq,<)i(-q,*)) ( A - 3 ) 

to Eq.( |A.l| ) the generating functional 

Wofij] = In JvhJ Vhexp {j [h, h] + Jj[h, h}} (A.4) 
is conveniently evaluated by solving the set of initial value problems given by 



2I)/;(q./J- (^ + i/q 2 )Mq./.i+./(q./)-<> : h(q.l) - 



d 

— -z^q 2 ) /i(q,t)+j(q,t) =0 ; fc(q, oo) = 



(A.5) 
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for h and h. The more general initial condition /i(q, 0) = /io(q) can be incorporated in the 
source field j(q, t) by the replacement 

J'(q,t)-»i(q,t) + *(*y»o(q)- ( A - 6 ) 

The solution of Eq.( |A.5| ), which is equivalent to calculating the inverse of the operator A 
(see Eq.( A.2[ )), is given by 



/i(q,t) J -A) V G (q,*,f) C (q,M') J \ j(q,t') 



where 



G (q,M') = eCt-Oe-^t*"*') and 



C (q,i,O = -^-f e -^ 2 l'-''l- e -^ 2 (*+*') 



i/q^ 

are the response and correlation functions of Gaussian theory for the KPZ equation, respec- 
tively. From Eqs.( [A.4| ) and ( [ASQ one obtains for the generating functional 

>v [j,i] = J J Q dt J q dt' (j(- q ,t)G (q,M )iM ) + 2i(-q>*) c o(q,M)j(q,* )J . 

(A.9) 

For the general initial condition /i(q, 0) = /io(q) the corresponding generating functional is 
obtained by applying the replacement Eq.( |A.6| ) directly to Eq.( |A.9| ). The response and 
correlation propagators can now be obtained by functional derivatives of Eq, (|A,8| ) with 
respect to j and j: 

fi?o(q, t\ q', = (/>(q, t)h(q[ , t')) = (27r) d <5(q + q')G (q, t, 

(A.10) 

Cb(q, t; q', t') = (/t(q, t)h(q' , t'))o = (2vr) d <5(q + q')0)(q, i, 0, 

where (. . .)o denote the average with respect to the Gaussian distribution generated by 
Eq.( A.l| ). From momentum conservation it is obvious that the full two-point correlation 



functions Q(q, t; q', t') and C(q, i;q',i') can be written in the same form as their gaussian 
counterparts (see Eq.( [A.10|) ) which serves as the definition of the full response function 
G(q,t,t') and the full correlation function C(q, t,t'). One should also note that the simul- 
taneous requirements /i(q, 0) = and h(q, oo) = forbid a Fourier transformation with 
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respect to time so that one has to stick to the above mixed representation of the propaga- 
tors for further calculations. Especially the normalization conditions imposed on correlation 



functions in order to define renormalized quantities (see Eqs.( |2.6| ) and (2.7)) have to be 
reformulated accordingly. Note that the exponents a and z implied by Eq. ([A.8 ) are the 
Edwards- Wilkinson exponents given by Eq.(|0|). 

In close analogy to the considerations described above the Gaussian part of the dynamic 
functional for Eq.( |1.9[ ) is the same as for model B of critical relaxation [17, 19] and can 
be written in the same symmetric form as given by Eq.(A.l) together with the conditions 
/i(q, 0) = and /i(q, oo) = 0. In this case the self-adjoint matrix operator A is given by 



A 



2D 



9 rA 



d_ 

dt 



(A.ll) 



where q = |q| is the modulus of the momentum vector q. The generating functional given 
by Eq.( [A.4[) is evaluated by solving the corresponding initial value problem for /t(q, t) and 
/i(q, t) (see Eq.( A.5[ )). The solution can be written in the same form as Eq.( [A.7| ), where 
instead of Eq.(|A.8|) one has 



G (q,t,0 = Q(t-t')e 



-viq 4 (t-t') 



C (q,t,O 



D 

v\q l 



and 



(A.12) 



for the response and the correlation function, respectively, of Gaussian theory for Eq.( |1.9| ). 
With Go and Co taken from Eq,( A,12| ) the corresponding response and correlation propaga- 
tors are again given by Eq. flA.10 ). We close this section by noting that the exponents a and 
z implied by Eq.( A.l^ ) are given by 



4 and a = (4 - d)/2 



(A.13) 



in contrast to Eq.( |l.ll ) 



Appendix B: Perturbation theory 

Due to the presence of strong coupling fixed points in the KPZ equation for d > 2 perturba- 
tion theory is only of limited value as compared to perturbation theory for model A critical 
dynamics, for example. However, some rigorous relations can be proved by analyzing the 
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c = 



q't' 



q,t 



(b) 




Fig. 2: (a) Graphical representation of the response propagator Q and the 



correlation propagator Cq (see Eq.( A. 10 )). (b) Graphical representation of the 



vertex of Eq.(|l.5|) (see Eq^^)). The wiggly lines represent the response field 
h(q,t) and the straight lines represent the deposition field h(q,t). 

building blocks of perturbation theory for response and correlation functions and therefore 
some details concerning perturbative calculations for Eqs.(1.5) and ( |1.9| ) will be described 
below. 



For the response and correlation propagators given by Eq.(AlO) we use the graphical 
representation shown in Fig.|2](a). The vertex and its analytical expression can be read off 
from Eq.( |2.2| ), they are shown in Fig.||(b). The momentum carried by the response field in 
Figj2|(b) is — qi — q2- Contributions to response, correlation, and vertex functions can be 
constructed from the elements in Fig according to the standard Feynman rules of dynamic 
perturbation theory || 17, [HJ. As a first example we analyze the response function G(q, t, t'). 
Any contribution to G from a perturbation expansion can be cast into the form of the block 
diagram shown in Fig]3|. According to the Feyman rules the first vertex contribution to an 
arbitrary diagram for G has to be arranged as shown in Fig.||. The remainder of the diagram, 
which is not neccessarily one-particle irreducible, is indicated by the shaded triangle and may 
be interpreted as an arbitrary contribution to the three-point vertex function. To lowest 
order this three-point vertex function is shown in Fig.0(b). From the explicit momentum 



20 




Fig. 3: Block diagram for the response function G(q, i,t'). The shaded 
triangle consists of an arbitrary number of vertices and propagators. To lowest 
order it is given by the vertex displayed in Fig.||(b) (see main text). 



(a) 




(b) 




Fig. 4: Block diagrams for the correlation function C(q, t, t'). (a) Incom- 
ing correlation propagator Co(qi, t\\ q', t'). (b) Incoming response propagator 



£/o(q')*';qi)£i) (see Eq.( |A.10 ) and main text). The shaded triangle has the 
same meaning as in Fig.|[ 

dependence of the vertex it is obvious, that for zero momentum q' = q the block diagram 
displayed in Fig.0 vanishes identically. One therefore finds the exact relation 



G(q = 0, t, t') = G (q = 0, t, t') = 9(t - t') 



(B.l) 



for the response function of Eq.(1.5). 

In contrast to Fig.|| the perturbative contributions to the correlation function C(q, t, t') 
cannot be represented by a single block diagram. Instead, two types of block diagrams are 
required as shown in Fig.||. Due to the initial condition /i(q, 0) = both block diagrams 
vanish identically for t' = 0. Following Ref.JT?]] Fig.|| is used to obtain an exact expression 
for the derivative of C with respect to the time argument t' . The diagrams for dC/dt' are of 
the same form as those for C. The main difference between the diagrams shown in Fig.||(a) 
and Fig.||(b) is that in (b) the internal time t\ is restricted to the interval < t\ < t' due 
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to causality, so that this block diagram vanishes identically for t' = 0. The remaining block 
diagram (a) is of the same type as the block diagram for the response function G shown 
in Fig.||. One therefore has a termwise correspondence between the perturbation series for 
dC(q,t,t')/dt'\ t/=0 and G(q,t,t' = 0). Gaussian theory (see Eq.flAlf)) yields 

^ C (q, t, t' = 0)= 2D e'^ 2t = 2D G (q, t, t' = 0) (B.2) 

and Figs|| and ||(a) then show that the two perturbation series only differ by an overall 
factor 2D. Therefore Eq.( [B.2D implies the relation 

^C(q,M' = 0)=2DG(q,M' = 0) (B.3) 

between the correlation function C and the reponse function G of the KPZ equation order 
by order in perturbation theory. According to the standard Feynman rules the arguments 
presented above for C and G also hold for arbitrary n-point correlation functions which 
differ only in the propagator (response or correlation) assigned to one of the external legs. 
Therefore Eq. flB.3j ) already establishes the proof of the operator identity Eq.( 2.1C| ) used in 
Sec. II. However, it must be pointed out here that the above arguments only constitute a 
rigorous proof of Eq^BJ?) and therefore of Eq.( [2.10| ) if the renomalization group fixed point 
is accessible by perturbation theory. For the KPZ equation this is only possible in d = 1. 
In d > 2 one encounters the well known strong coupling behavior which forms a formidable 
obstacle for analytic theories of dynamic scaling of Eq. ( |1.5|) || |l0| . For the above derivation 
this means that Eq.( [B.3D may not hold in d > 2 at the renormalization group fixed point 
despite its validity to all orders in perturbation theory. Nonetheless the above perturbative 
analysis provides some evidence that Eq.( B.3 ) and therefore Eq.( |2.10| ) holds beyond d = 1. 

For Eq. (|1.9| ) the building blocks of the perturbation theory can again be taken from Fig.||, 
with the modification that the response and the correlation propagator (see Eq.( A.lCj )) are 
now given by Eq.( A.12] ) and that the expression Aiqi -q2(qi +q2) 2 must be assigned to each 
vertex as can be read off from Eq.( |3.2| ). It is then straightforward to see that the arguments 
given above for the KPZ equation can be directly applied to ideal MBE dynamics, where no 
strong coupling behavior is encountered in any spatial dimension of physical interest. The 
exact relations given by Eqs.(B.l) and ( |B.3[ ) and the operator identity Eq.(2.10) therefore 
also hold for Eq.O). 
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Appendix C: Response and correlation functions 



In order to justify Eqs. Q2.7p and Q3.6| ) within the dimensional regularization scheme || in the 
q, t representation and to obtain some indication how the short-time to long-time crossover 
takes place the response and correlation functions of Eqs.(L5) and Ql.9|) are calculated here 
to one-loop order. 

The one-loop contribution to the response function for the KPZ equation is given by the 
block diagram shown in Fig|], where the shaded triangle is replaced by a single vertex shown 
in Fig.|2|. The analytic expression for this diagram is then given by 

d d q' 



Gi(q,t,0 



dU 



dto 



(q • q)(q • (q - q)) 



10 JO J (27T) d 

x G (q, t u O^oCq' - q, *2, h)C (q', t 2 , ti)G (q, t, t 2 ), 



(C.l) 



where Go and Co are given by Eq.(A.8). For simplicity we only consider Eq.( p.ip in the 
limit q — > so that we can employ the expansion 



G (q! - q, t 2 ,h) = G (q, t 2 ,h) 1 + 2i/(q' • q)(t 2 - h) + 0(q 



(C.2) 



The q' integration in Eq.( C.l ) to leading order in q is then reduced to the calculation of 
second moments of a Gaussian in d dimensions. The result is 



Gl(q,*,t') 



q 



t rt 2 

dt 2 / dt\ 
t' Jt' 



G^t,t'){2u) 2 - d l 2 
d-2 



(C.3) 



d/2 



d-2 



h \ .-d/2 

6o 



2d y ' V 2d 2t 2 , 

where the effective coupling constant g is defined by Eq.( [2.5p . The remaining integrals in 
Eq.(C.3) can be easily performed using dimensional regularization Q with d = 2 + e in the 
exponents of t 2 —t\ and t 2 . Note, that the prefactor d — 2 in Eq.( p.3| ) comes from an angular 
integration and must not be cancelled by factors 1/e indicating UV singularities in the time 
integral Q. With the definition of u according to Eq.(^j]) (Z g = 1 at this order) one obtains 
for G(q, t, f) = G (q, t, f) + Gi(q, t, f) 



G(q,t,t') =G (q,t,t')l 1- \uijr E 



d ~ 2 (o„n _ ,'^-d/2 ^d-A (2u(t - t')f 
de { [t t)] + 4d (2vt) d / 2 



(C.4) 

The 1/e pole (the UV singularity) in Eq.( C.4 ) can be removed, e.g., by requiring G(q, t R ,0) 
to stay finite for e — > 0, where t R = 1/ (2p?v R ) is a reference time and u R is given by 
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Eq. Q2.6| ). Minimal subtraction yields the renormalization factor Z v quoted in Eq.( |2,7| ). Note 
that the short-time contribution to G does not produce an additional 1/e pole. By naively 
exponentiating the q-dependence of G in the long-time limit one obtains at the infrared 
stable fixed point u = u* ^ 

.2- 



G R (q,t,t') = 0(t-t')exp 



(C.5) 

The predictive value of Eq.( |C.5| ) is very limited because u* is infinite for d > 2. In d = 1 
Eq.( C.5 ) indicates that the combination (i — t') 2 jt d l 2 of the time arguments governs the 



crossover to the long-time scaling behavior of G for t — » oo with fixed t—t'. From dimensional 
arguments and the fact that the short-time contribution to G does not produce additional 
UV singularities we can infer that according to Eq.( |C,5| ) q 2 (i — t') 2 jt d l 2 for d = 1 is the 
scaling argument which governs the leading finite-time correction to the asymptotic long- 
time behavior of G. Furthermore, Eq.( |C.5| ) shows that G R is analytic in t' for t' <C t at 
the one-loop level, but this behavior may be modified in higher orders. Finally, we note 
that the scaling form of the asymptotic long-time contribution to G R (q,t,t') given by the 
exponential in Eq.( |C.5| ) has recently been derived by combining perturbative methods with 
a mode coupling theory for the KPZ equation jin]. 

The correlation function C(q, t,t') for Eq.( |1.5| ) can be discussed in much the same way 
as the response function. This time we simplify the calculations even further by limiting 
ourselves to q = 0. In this case only the diagram in Fig.|](b) contributes and we obtain to 
one-loop order 

C(0, t, = 2D min(t, t') + y jT dt 2 J* dt x J -^q A [C (q, h , t 2 )} 2 , (C.6) 



where Cq is given by Eq.( A.8 ). The integrations in Eq. (|C.6|) can be easily performed and 



using dimensional regularization one arrives at 

C(0,t,t'<t) = 2£>jt' + ^i^ [{2v{t-t')) 2 - d l 2 - {2v{t + t')f~ d/2 (C.7) 
+ {2vtf- d l 2 ((4 - d)(t'/t) - d{t' /t) 2 "^ 2 )] } , 

where Eq.fl2.6l) has been used with Z g = 1. The 1/e pole in Eq. flC.7D can be removed by 
demanding C(0,t R ,t R ) = finite, where t R = l/(4:fj, 2 u R ) is chosen as the reference time. 
Using minimal subtraction one finds the renormalization factor Zd quoted in Eq. (12.71). For 
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if <C t Eq.( p,7| ) can be simplified to 

C(0,t,t' < t) = 2Dt' 



i_!f^( 2 ^V) 1 - d / 2 + o(t' 2 ; 

S 2 



(C.8) 



which explicitly shows that the short-time contribution to G produces an additional 1/e 
pole. For t' > the renormalized correlation function can be naively exponentiated at the 
infrared stable fixed point u = u* ^ 0. The result is 

C R {0,t,t'<t) = D R (t + tf) (2^V(i + i'))^ 1 " (*-0 (2^V(*-0) ei (C9) 



2et'(2u R fi 2 t) e 1 + 2t'(2u R fJ 2 t' 



where is the short-time exponent given by Eq.(2.18) and d = 1 has been assumed. The 
short-time scaling behavior for tf <C t is also reproduced by Eq.(C.9). However, from 
Eqs.( p.lOj ) and (B.l) one expects ^C R (0,t,t' = 0) = 2D R which is not reproduced by 



Eq. (|C.9[) , because 9 > 1. Therefore Eq. (p.9|) can only give a rough idea of the true scaling 



form of the correlation function C R for the KPZ equation. However, Eq.( C.9 ) indicates, 
that for q = short-time corrections to the correlation function persist indefinitly (see also 

Eq.dazD). 

For ideal MBE dynamics according to Eq.([0]) the one-loop contribution to the response 
function is again given by the block diagram shown in Fig.^, where the shaded triangle is 
replaced by a single vertex. The analytic expression for this diagram is then given by 

Gi(q,t,0 = 4A 2 q 2 y o dt x J Q dt 2 J S-(q'.q)(q'.(q'-q))(q'- q )2 ( C .10) 
x G (q, t x , t')G (q! - q, t 2 , h)C {q', t 2 , *i)G (q, t, t 2 ), 

where Go and Go are given by Eq.( A.12 ). For simplicity we only consider Eq.( C.10| ) in the 
limit q — > 0, i.e., we use the expansion 



G (q' " q. *2, tl) = G (q', t 2 ,h) [l + 4i/ iq ' 2 (q' • q)(t 2 - tx) + C(q 2 ; 
The q' integration in Eq.( |C.10|) to leading order in q yields 

GM = ^^ (q,M')(2^ 



(C.ll) 



(C.12) 



dU 



t2 



dU 



d-6 
d 



d-6 h 



-d/4 



where the effective coupling constant g\ is defined by Eq.(3.4). As in Eq.( |C.3| ) the remaining 
integrals in Eq,( C.12| ) can be performed using dimensional regularization with d = 4 + e 
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in the exponents of ti — t\ and t%. As usual the 1/e poles indicate UV singularities in the 
time integral. With the definition of u according to Eq.( |3.5| ) and Z gi = 1 one obtains for 
G(q, t, t') = Go(q, t, t') + Gi(q, t, t') in the limit t — > oo with t — t' = const 



G(q,t,0 = Go(q,*,0 l-y¥" E 



de " 4 4d (2u 1 t) d / 4 

(C.13) 

up to terms 0((t - t') 3 /t d ^ +1 ). The 1/e pole (the UV singularity) in Eq. flCl3D can be 
removed by the minimal subtraction scheme described above, where t R = l/(2/i 4 z/^) defines 
the reference time and vf 1 is given by Eq.(p,5|). One obtains the renormalization factor Z Vl 
quoted in Eq.( |3.6D . By naively exponentiating the q-dependence of G in the long-time limit 
one obtains at the infrared stable fixed point u = u* ^ 

.2" 



(™(q,i,t') = 9(t-t')exp 



3 g 4 d _ 8 ( 2 ^V(t 
1 — -u 



to,,",,*,/ - " i__„ 

4 ^ 8 d (2v^H) d/A 

(C.14) 

In contrast to Eq.([1.5|) the infrared stable fixed point for Eq.(|l.S|) is finite in any dimension 
of physical interest. Especially one has u* = 0(e) so that an e-expansion around the upper 
critical dimension d c = 4 can be performed. The qualitative behavior of G R according to 



Eq.( C.14| ) is very similar to the behavior of G R for the KPZ equation in d = 1 (see Eq.( p.5|) ). 



Here the leading finite-time correction to the asymptotic long-time behavior is governed by 
the combination (t — t') 2 /t d / i of time arguments. 

The one- loop contribution to the correlation function C(q = 0,t, if) for Eq.( |1.9|) vanishes 
identically due to an additional factor q 2 in the vertex (see Eq.( C.10| )) so that 



C R (0, t, t') = 2D min(t, t') + 0{u 2 ). (C.15) 



Eq.( C.15| ) directly demonstrates that Z& = 1 as quoted in Eq.(p^) to one-loop order. As 



in the case of KPZ dynamics Eq. QC.15 ) demonstrates that finite-time corrections to the 



correlation function for Eq.(|1.9|) persist indefinitely for q = 0. 
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